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SYSTEM AND METHOD FOR FAST CODE PHASE 
AND CARRIER FREQUENCY ACQUISITION IN GPS RECEIVER 

This application claims the benefit of U.S. Provisional AppHcation No. 60/223,973, filed 
August 9, 2000. 

BACKGROUND 

Field of the Invention 

[0001] The present invention relates generally to GPS positioning systems. More specifically, 

the present invention relates to improved code phase and carrier acquisition in GPS positioning 
systems. 

=0 Background of the Invention 
[OOOip The global positioning system (GPS) is employed to determine position in a multitude of 

g applications. For example, in navigation applications, a user uses a GPS receiver to determine 
. her instantaneous position as well as her position over time. Another use for GPS receivers 
I D tracking objects. For example, placing a GPS receiver in a truck allows tracking of the tmck's 
movements. Attaching the receiver to cargo allows precise tracking of the locations of that 
cargo. These concepts can be extended to fleet management for common carriers and delivery 
management for delivery service providers. Numerous additional and new services are being 
deployed and developed for using the position-determining capabilities of GPS receivers. 
[0003] Figure 1 is a schematic diagram of a conventional GPS receiver 101 used for determining 

position. In Figure 1, GPS receiver 101 is simplified to point out the primary relevant functions 
of a conventional GPS receiver. An antenna 102 receives a GPS signal fi-om GPS satellites 103a, 
103b, 103c and 103d. Antenna 102 applies the received GPS signal to a signal conditioning 
processor 104. Signal conditioning processor 104 amplifies, filters and downconverts the signal 
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to baseband for processing. The baseband signal is applied to carrier and code phase tracking 
algorithms in processing block 106. Processing block 106 contains a multiplier 108, a correlator 
1 10, a carrier frequency oscillator 1 12, a Gold code generator 1 14 and an integrator 1 1 6. 
Multiplier 108 multiplies the baseband signal by an estimated carrier frequency received from 
carrier frequency oscillator 112. Carrier frequency oscillator 112 can be a voltage controlled 
oscillator (VCO) or a numerically controlled oscillator (NCO). Correlator 110 correlates the 
signal with a replica of a Gold code generated by code generator 1 14. The Gold code is a unique 
and known code generated by each GPS satellite. The terms "code" and "Gold code" are used 
interchangeably herein. The output of correlator 110 is integrated in integrator 1 16. The output 
1% of integrator 1 16 is input to a digital signal processor 1 18 to generate information required for 

[f code tracking generator 1 14 and carrier frequency oscillator 112. This information includes 

ill 

,g carrier phase and code phase information. 
[0004]^ The combination of carrier frequency oscillator 1 12, multiplier 108, integrator 116 and 

^ Jf DSP 118 operates as a carrier tracking loop that compensates for errors in the down converter 
]^ frequency reference as well as any Doppler shift associated with a particular satellite. The 
combination of code generator 1 14, correlator 110, integrator 1 16 and DSP 1 18 operates as a 
code tracking loop that compensates for any shift in the code phase of the received signal. 
[0005] Generally, there are a series of processing blocks 1 06 running in parallel in a 

conventional GPS receiver. Each of the series of processing blocks is tuned to a different GPS 
satellite. Generally, there are at least 2 processing blocks 106 per GPS satellite. 
[0006] GPS receiver 101 performs operates in two modes to process a GPS signal. First, the 

receiver must acquire the GPS signal's carrier frequency and Gold code phase. This is known as 
the acquisition mode. Second, receiver 101 must track the carrier frequency and Gold code 
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phase using the carrier frequency and Gold code phase determined in the acquisition phase as a 
starting point. This is known as the tracking mode. 
[0007] During the acquisition mode, a GPS receiver uses a priori knowledge. For example, the 

receiver knows the nominal carrier frequency of the GPS signal. However, the nominal carrier 
frequency is not likely to be the one actually received by the receiver for several reasons. First, 
Doppler shifl caused by the relative motion of the transmitting satellite and the receiver changes 
the received carrier frequency. Moreover, the local oscillator (not shown) introduces its own 
en-ors. Consequently, both the carrier frequency and offset of the Gold code (Gold code phase) 
must be determined. 

[O0O8| To make these determinations, conventional GPS receivers generally perform a serial 

'5 search by trying different combinations of carrier frequency and code phase until the output of 
integrator 1 16 exceeds a pre-determined threshold, indicating that a match has been obtained. 

\ y 

3 The search must be performed for each GPS satellite the receiver will use to determine position. 
?fl Using multiple processing channels as described above, the GPS receiver can perform the 
searches in parallel. 

[0009] Figure 2 is a graphical representation of a search process performed by conventional GPS 

receiver 101 to determine carrier frequency and Gold code phase. Oscillator 1 12 is typically set 
to the nominal carrier frequency of the GPS signal. In addition, the code is generated by code 
generator 1 1 4 with an initial code phase. If the received carrier frequency and code phase match 
the carrier frequency generated by oscillator 112 and the code phase of the code generated by 
code generator 114 respectively, the output of integrator 1 16 is relatively high. On the other 
hand, if the carrier frequency or the code phase do not match, the output of integrator 1 16 is 
relatively low, and may be essentially noise. 
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[0010] Generally, a conventional search algorithm first sets a frequency and lets the Gold code 

phase vary. As the code phase varies, controller 1 1 8 measures the output of integrator 116. If 
the output of integrator 116 exceeds a pre-determined threshold, controller 118 presumes that 
lock has been achieved and controller 118 places the receiver in its tracking mode using the 
values determined by the search algorithm for carrier frequency and code phase. If the output of 
integrator 116 does not exceed the pre-determined threshold, controller 1 18 presumes that lock 
has not been achieved. Once all code phases (code phase cpl through cpN) have been tried, 
controller 118 switches carrier frequency oscillator to output a carrier fi-equency f2. The process 
is repeated for each code phase cpl through cpN, and each frequency frequency 1 through 

^0 frequency M, unless controller 1 1 8 stops the process early because the output of integrator 1 1 6 

bp 

indicates a lock condition. If no frequency code phase combination is found, the receiver cannot 
acquire the signal. 

[001 if Conventional GPS receivers perform each step in the above process on a different part of 

m the received signal. That is, at each step, a new portion of the GPS signal received by the 

:y receiver is processed. While conventional GPS receivers can sometimes acquire the GPS signal 

P 

H in a few seconds, in general acquisition is significantly longer. 
[0012] After the acquisition mode, the receiver switches to a tracking mode. In the tracking 

mode, the receiver continually adjusts the frequency of oscillator 112 and the offset of the Gold 
code phase to maintain lock to the incoming carrier and code phase. Code phase tracking is 
typically assisted in a well-known maimer using early and late code generators that generate 
early and late codes respectively (as well as a punctual code). If the received code correlates 
better with the early code, controller 1 1 8 reduces the phase delay of code generator 114. If the 
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received code correlates better with the late code, code generator increases the phase of code 
generator 114. No change is made if the received code correlates best with the punctual code. 
[0013] Carrier tracking is typically accomplished by analyzing the phase of the signal output by 

the integrator. If the carrier frequency generated by oscillator 1 12 matches the received carrier 
frequency, then there will be no phase rotation. If there is no such match, then controller 118 
increases the frequency or decreases the frequency of oscillator 1 18 to compensate for any phase 
rotation. 

[0014] In addition to tracking, the receiver can demodulate the 50Hz data modulated on the 

carrier of the GPS signal. This 50Hz signal carries information related to the particular satellite 
that transmitted the signal. For example, this signal contains ephemeris information that 
provides information on the position and motion of the GPS satellite. 



in 



SUMMARY OF THE INVENTION 

The present invention is related to improving the acquisition phase of a GPS receiver by 
reducing the time required to acquire the GPS signal, using only a very short portion of the 
received GPS signal. That is, the present invention reduces the time required for the receiver to 
determine the proper carrier frequency and Gold code phase. This information can replace the 
acquisition phase conventional GPS receiver by providing the initial value of carrier frequency 
and Gold code phase for the receiver to use in its tracking phase. In addition, the code phase 
measurements derived using the present invention can be used themselves to determine position 
in a system such as the system described in U.S. Patent No. 6,243,648 to Kilfeather et al,, which 
is hereby incorporated by reference herein in its entirety. 
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[0016] In one embodiment of the present invention, a GPS signal is received and conditioned by 

a signal-conditioning processor. The GPS signal so received is a composite signal having 
contributions &om all of the GPS signals in view of the receiver. A short portion of the signal is 
collected and stored in a memory. An FFT is performed on the short portion of the signal to 
convert it to the frequency domain. A Gold code sequence is generated. An FFT is performed in 
the generated Gold code sequence to convert it to the frequency domain. Because the Gold code 
is a known sequence, generating the Gold code sequence and converting it to the frequency 
domain are preferably pre-operability events. That is, the Gold code sequence and its frequency 
^ domain representation can be pre-computed prior to using the GPS receiver to determine 

y 

% position. In that manner, the Gold sequence and/or its frequency representation are stored in the 

iiJ 

p memory for ftiture operation. The FFT of the short portion of the input signal and the Gold 

in 

,p sequence are multiplied. An inverse FFT operation is performed to convert the product to the 

lU 

time domain. The time domain representation is the convolution of the short portion of the input 
signal with the generated Gold code sequence. A magnitude calculator calculates a point-by- 
;=? point magnitude of the time domain product. A peak detector determines where the peak is 
located. The peak represents the code phase in the short portion of the signal. 
[0017] If no peak is present, the controller moves to a new carrier frequency and repeats the 

process. This can be performed by multiplying a time domain signal by a complex exponential 
having a frequency equal to the frequency shift. Alternatively, this can be performed by shifting 
all bins in the frequency representation of the input GPS signal or Gold code one way or the 
other by a predetermined number of bins. The frequency shift will be equivalent to the 
resolution of the FFT times the number of bins shifted. The Gold code detection process is then 
continued. This continues until a peak is observed, or if no peak is observed after all desired 
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frequencies have been tested, the process ends without successful acquisition. The entire process 
is repeated for each GPS satellite contributing to the input signal by using the Gold code unique 
to each such GPS satellite. 
[0018] In a another embodiment of the present invention, the process described above is 

performed on a longer duration of the input GPS signal so that integration can be used to make 
the code phase peak more distinct. In this embodiment, a multiple millisecond portion of the 
input GPS signal is collected and stored in memory. The input GPS signal is divided into one 
millisecond segments. Each millisecond segment is converted to the frequency domain as a 
frequency domain representation of that segment. The frequency domain representations are 
^0 multiplied by a frequency representation of the Gold code unique to one of the satellites in view 
of the receiver. Preferably, the frequency representation of the Gold code was generated and 
stored prior to operation of the receiver. The product of each multiplication is converted to the 
[ time domain and represents the convolution of the each of the 1 ms segments with the Gold code 
m for the particular satellite. The magnitude of each point in the convolutions is calculated. The 
^0 magnitudes are integrated (summed) point-by-point across the convolutions. A peak is then 
H detected. The location of the peak corresponds to the code phase. The carrier frequency is 

determined as above by appropriate shifting of the input signal or the Gold code in the frequency 
domain. 

[0019] In another embodiment of the present invention, the integration discussed above is 

performed using coherent processing. Two steps are performed. The first step is to determine a 
coarse code phase (the peak location in the convolutions) and a coarse carrier frequency using 
the multiple millisecond method described above. Then, a complex sine wave is fitted to the 
points at the location of the peaks in the convolutions. The frequency of this sine wave is the 
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difference between the received carrier and the estimated coarse carrier frequency. The points in 
the convolutions are phase rotated to remove this frequency difference by complex multiplication 
of each convolution with a complex exponential at a frequency equal to the fitted sine wave. 
Doing so makes the convolutions coherent to one another. The convolutions are then integrated 
(summed) point-by-point. The magnitude of the sum is taken at all of the points and the location 
of the peak is detected. The process is optimized by rotating, summing and taking the magnitude 
of only a few points around the probably location of the peak, as well as the peak itself. 
Interpolation techniques are then applied to further refine the location of the peak. 
[0020] The present invention is described in greater detail in the detailed description of the 

invention, the appended drawings and the attached claims. 

■| BRIEF DESCRIPTION OF THE DRAWINGS 
[0021] Figure 1 is a schematic diagram of a conventional GPS receiver used for determining 

IB position. 

5 ^ 

[0022p Figure 2 is a graphical representation of a search process performed by conventional GPS 

P 

i 

receiver to determine carrier frequency and Gold code phase. 
[0023] Figure 3 is a schematic diagram of a GPS receiver for acquiring a GPS signal received 

from one or a plurality of GPS satellites that can be used according to the present invention. 
[0024] Figure 4 is a schematic diagram of a processor for determining the location of the Gold 

code sequence. 

[0025] Figure 5 A is a schematic diagram of frequency domain processing according to a first 

embodiment of the present invention. 
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[0026] Figure 5B is a schematic diagram of frequency domain processing according to a second 

embodiment of the present invention. 
[0027] Figure 6 is a schematic diagram of an embodiment of a system for processing a multiple 

millisecond portion of the input GPS signal according to a third embodiment of the present 

invention. 

[0028] Figure 7 illustrates coherent processing in accordance with the present invention. 

[0029] Figure 8 illustrates the use of a table that stores pre-computed timing offsets of the input 

signal relative to the Gold code. 

DETAILED DESCRIPTION OF THE INVENTION 
[0030p The present invention is related to improving the acquisition phase of a GPS receiver by 

reducing the time required to acquire the GPS signal, using only a very short portion of the 
'JJ received GPS signal. Preferably, the short portion is between a millisecond and a second of the 
I " input GPS signal. Thus, the present invention reduces the time required for the receiver to 
m determine the proper carrier frequency and Gold code phase. The invention allows for signal 

£ ""1 

acquisition using collected portions of the signal as short as 1 millisecond in duration. 
[003^]"'= Figure 3 is a schematic diagram of a GPS receiver 301 for acquiring a GPS signal 

received from a one or a plurality of GPS satellites that can be used according to the present 
invention. For the sake of clarity only a single charmel will be described in the following. 
However, it would be apparent to those skilled in the art that the single channel described can be 
replicated with suitable changes in Gold code generation to process the input signal for all GPS 
satellites in view simultaneously. Each channel uses a Gold code generator to generate a Gold 
code corresponding to the particular GPS satellite whose signal the charmel is intended to 
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receive. Alternatively, a single channel is used in a serial mode by processing the signal using 
the Gold code for each GPS satellite in turn. 
[0032] Receiver 301 includes an antenna 302 that receives a composite GPS signal That is, the 

input GPS signal received by receiver 301 has contributions from all of the GPS satellites in 
view of the receiver 301 . Antenna 302 is coupled to a signal conditioning processor 304. Signal 
conditioning processor 304 amplifies and filters the GPS signal collected by antenna 302. 
Preferably, signal conditioning processor 304 also converts the signal to some other frequency, 
such as baseband, for processing. The baseband signal is digitized by an analog-to-digital A/D 
converter 306. The digitized signal is then stored in a memory 312. Memory 312 is accessible 
=0 to, or part of, a digital signal processor (DSP) 308. A processor 310 using memory 312 

processes the GPS signal as described below to acquire the signal, Le., to determine the carrier 

rat 

' Ji frequency and code phase for each of the GPS satellites contributing to the GPS signal. 
[0033]" A/D converter 306 is any A/D converter that can digitize the signal at twice the 

m information bandwidth. In this case, the information bandwidth is the GPS chip rate of 1 .023 

MHz (assuming signal-conditioning processor 304 downconverts the signal to baseband). For a 
\^ complex signal, therefore, the minimum sampling rate is 1 .023 MHz. In general, a higher 

sample rate is better, but the sample rate is constrained by the size of memory 312. Preferably, 

the sample rate is 2.048 MHz. 
[0034] Any portion of the GPS signal can be stored in memory 3 12. Preferably, however, the 

portion of the GPS signal stored is more than millisecond and less than a second. More signal 

allows longer integration times, which is useftil where the signal-to-noise ratio (SNR) is low, but 

comes at the expense of a longer processing time. This is a design trade-off that is determined 

by each implementation of the present invention. 
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[0035] As mentioned above, the digitized GPS signal contains contributions from each of the 

GPS sateUites in view of the receiver. The receiver determines which of the satelHtes is in view 
using pubHcly available information in a number of well-known ways. For example, the receiver 
can correlate the received signal with all published GPS Gold codes. Alternatively, if the 
receiver has a current almanac, the current time and approximate location (e.g., from its last 
position determination), the receiver can compute the GPS satellites in view. Further, the GPS 
receiver can be told which satellites are in view from a communication link. In this manner, the 
receiver knows which Gold code sequences to apply in its processing. 

[0036] The Gold code phase and carrier frequency are determined from the stored input GPS 

'0 signal as follows. The terms "delay of the Gold code," "Gold code phase" and "code phase" are 
'z, used interchangeably herein. As mentioned above each GPS satellite is assigned a unique Gold 
code that is transmitted in the GPS signal. Thus, the composite signal collected by antenna 302 
, contains the GPS signal, and the corresponding Gold code for each satellite in view of the GPS 
in receiver. The Gold codes are publicly available. 

[0037p There are 1 ,023 bits or chips in the Gold code, which are transmitted at a rate of 1 .023 

Mhz. Thus, a complete Gold code sequence has a 1-ms duration (i.e., repeats every millisecond). 
Every 20 Gold code sequences (20 milliseconds) there may be a phase shift (or transition) of the 
GPS signal by 180 degrees. If a phase shift occurs, the signal inverts. Consequently, obtaining 
at least a 1-ms portion of the GPS signal guarantees that at least one fall sequence of the Gold 
code is captured. The problem then is to find where the sequence begins in the 1 millisecond 
segment of the input GPS signal that has been collected. 

[0038] Figure 4 is a more detailed schematic diagram of processor 3 10 for determining the 

location of the Gold code sequence. The digitized signal stored in memory 312 is convolved 
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with the known Gold code to determine where the Gold code begins. The known Gold code is 
generated by a Gold code generator 402. The Gold code generated by Gold code generator 402 
is convolved with the received GPS signal in convolver 404. Convolver 404 outputs a convolved 
signal 406. Each peak 406a, 406b, 406c or 406d represents where the Gold code generated by 
Gold code generator 402 matched best with the Gold code in the digitized GPS signal. Note that 
peak 406c is inverted. This is an example of the effect of the 180 degree phase shift (also 
referred to as "transition" herein) inherent in the GPS signal that occurs potentially every 20 
milliseconds. 

[0039]_ As is well-known, convolution in the time domain is a computationally intensive process. 



Consequently, the present invention preferably uses frequency domain techniques for performing 



, embodiment of the present invention. For clarity, the processing is described for a 1-ms portion 
iQ of the input GPS signal. Gold code sequence generator 402 generates a 1-ms (Le., complete 

5 I 

cycle) portion of the Gold code. The Gold code is input into an FFT algorithm 506a. The 

U 

digitized 1-ms portion of the received GPS signal 504 is input to an FFT algorithm 506b. The 
output of the FFT algorithms is the frequency domain representation of the signal. The 
frequency domain representations of the input GPS signal segment and the Gold code are 
multiplied by a frequency domain multiplier 508. 
[0041] As is well-known, multiplication in the frequency domain is equivalent to convolution in 



the convolution of the received GPS signal with the known replica of the Gold code. 




Figure 5A is a schematic diagram of frequency domain processing according to a first 



the time domain, but requires far fewer operations. The frequency domain multiplication 



product is converted back to the time domain by an inverse FFT algorithm 510. A magnitude 



calculator 512 performs a point-by-point magnitude calculation on the output of inverse FFT 
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algorithm 510. A peak detector 514 determines the location of the peak by determining where a 
bin value exceeds a pre-deteraiined threshold. That bin location is taken as the code phase of the 
received GPS signal. Interpolation techniques can be used to refine the peak location and are 
described below. 

[0042] Any size FFT can be used. Hov^ever, in one embodiment of the present invention, the 

FFT resolution varies depending on which signal is being processed and what stage the process 
is in. Moreover, FFTs are not required and powers of 2 are not required. However, powers of 2 
are preferred as they permit use of computationally efficient algorithms such as the FFT. In one 
embodiment of the present invention, the FFT resolution is 1 KHz/bin. The input signal is 
digitized at 4.096 MHz. To obtain the IKHz resolution a 4,096 point FFT is used to convert the 
incoming signal to the frequency domain. Because the information bandwidth is 1.023 MHz, the 

!C 

signal is oversampled. Consequently, only the lowest 1 ,024 points are used, yielding a 1 .024 
, MHz effective bandwidth FFT with a 1 KHz resolution. 
[00435 As described above, the Gold code sequence is unique for each GPS satellite. It is a 

^0 permanent and unchanging code. Consequently, the FFT of the Gold code sequence can be pre - 

O 

computed and stored in a non-volatile memory, such as ROM, PROM, EPROM, EEPROM or 
any other non- volatile memory. In addition, the FFT can be stored in memory 312 upon 
initialization. Figure 5B is a schematic diagram of a system according to a second embodiment 
of the present invention that can be used when the frequency domain representation of the Gold 
code is pre-computed and stored prior to operation of the system. The pre-computed frequency 
domain representation of the Gold code 520 is input to frequency domain multiplier 508 and 
processing continues as described above. 
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Another improvement offered by the present invention arises from recognizing that the 
contribution of each GPS satellite in view of the GPS receiver is present in the composite input 
GPS signal. Consequently, the same signal portion can be used for all GPS satellites. 
Consequently, the input signal segment needs to be converted to the frequency domain only one 
time. 

As mentioned above, if the carrier frequency used to mix the signal to baseband is not the 
true carrier frequency, the output of multiplier 404 or 508 is essentially noise. Consequently, 
there will be no discemible peak in the magnitude generated by magnitude calculator 512. That 
is, there probably will be no bin that has a value greater than the pre-determined threshold for the 
peak. When this occurs, processor 310 assumes that the carrier frequency used to generate the 
Gold code is incorrect. Thus, it changes the carrier frequency by a pre-determined amount and 
repeats the above process. It continues repeating the above process until a set of discrete carrier 
frequencies is tested that covers the expected Doppler uncertainty and the frequency uncertainty 
of oscillator used to mix the signal to baseband. Once an carrier frequency is found that causes a 
peak (a bin having a value greater than the pre-determined threshold), processor 310 discontinues 
searching for a carrier frequency and uses the carrier frequency that provides a peak. 

Shifting the frequency using the present invention is relatively simple since the digitized 
signal has already been converted to the frequency domain. This is because the frequency of the 
signal is changed simply by shifting the FFT of the signal by N bins. The total frequency shift is 
N times the resolution of each bin. Shifting to the left decreases frequency, while right shifting 
increases frequency. If a frequency shift of other than the bin resolution is desired, it is 
accomplished by multiplying the time domain representation of the signal by a complex 
exponential having a frequency equal to the desired offset prior to performing the FFT. 
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[0047] In addition, to minimize losses due to not matching the actual carrier frequency, the 

carrier frequency is also stepped in half-bin resolutions. The way this is accomplished is to 
multiply the digitized signal stored in memory 312 by a complex exponential having a frequency 
equal to half a bin resolution. The FFT of this half-bin shifted frequency is then taken and the 
above procedure is repeated for this signal. The highest resulting peaks are compared, and the 
carrier frequency and code phase corresponding to the higher of the two peaks is chosen as the 
code phase and carrier frequency. It should be noted that better results are obtained for smaller 
frequency shifts, Le.^ higher resolution. However, such higher resolution processing requires 
more processing time to perform. 

[004Sj The Gold code is sampled at 1 .024 MHz. Although the Gold code could be sampled at 

its inherent 1.023MHz, it is preferable to sample at a rate that is a power of 2, so that efficient 
algorithms such as the FFT can be used. A 1024-point FFT is applied to the sampled Gold code 

llj 

. sequence. 

[0049p As mentioned above, the Gold code sequence is unique to each satellite. It is permanent 

{i.e., does not change for a given satellite) and well-known. Therefore, the Gold code and 
consequent frequency representation for each satellite can be pre-computed and stored in a table 
for later access. Thus, the Gold code sequences for each satellite are applied to FFT algorithms 
to generate frequency representations of the respective Gold codes. These frequency 
representations of the Gold codes are stored in a memory before the receiver is placed in 
operation. The frequency representations can be stored in any type of non- volatile memory, for 
example ROM, EPROM, EEPROM and any other memory device accessible by processor 310. 
Moreover, the frequency representations are preferably stored in a table. When processor 310 
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requires the frequency representation of a particular satellite's Gold code, it retrieves it from the 
table, rather than computing it at runtime. 

The pre-computed Gold code sequence also provides more efficient performance of the 
half-bin analysis described above. In this case, the frequency representation of the Gold code is 
shifted by a half bin and stored in memory as described above for the non-shifted Gold code. 
This shifted Gold code is used for the half bin analysis. The improved computational efficiency 
must be weighed against the additional memory resources required to store the pre-computed 
half-bin shifted Gold code. As with shifting the input signal in frequency, better results are 
obtained for smaller frequency shifts, /.e., higher resolution. Again, however, such higher 
resolution processing requires more processing time to perform. 

The frequency domain representations of the input signal and the Gold code for a 
particular satelUte are multiplied against one another to generate a product signal. The product 
signal is the frequency domain representation of the convolution of the input signal and Gold 
code. Inverse FFT algorithm 510 performs the inverse FFT operation to obtain the time domain 
convolution of the two signals. 

Inverse FFT algorithm 5 1 0 could perform a 1 024 point inverse FFT. However, the 
results can be substantially improved by increasing the resolution of the inverse FFT operation. 
This is done by appending O's to the end of the product signal. For example, doubling the 
resolution is achieved by appending 1024 O's to the end of the product signal and taking a 2048- 
point inverse FFT. Likewise, quadrupling the resolution is achieved by adding 3072 O's to the 
end of the product signal and taking a 4096-point FFT. Increasing resolution of FFTs and 
inverse FFTs in this manner is well-known to those skilled in the art. 
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[0053] Increasing the resolution of the inverse FFT operation reduces the risk of the real peak of 

the convolution of the input signal and Gold code falling between bins of the inverse FFT 
operation. Spreading the peak energy between bins can make the peak harder to detect, 
especially in the high noise environment in which GPS receivers operate. Moreover, the peak 
energy in either bin may not be greater than the pre-determined threshold. As a result, peak 
detector 514 will not detect a peak that is validly there. Increasing the resolution of inverse FFT 
operation 510 avoids this problem. 

[0054] Another consequence of taking a 4,096 point FFT of the input signal is that when 

performing the frequency shift described above, the shifting should be performed using all 4,096 
points in the original FFT. Although the shift could be done using just the 1,024 lowest points, 
more accurate results are obtained by shifting using the 4,096 point FFT, and using the lowest 

i P 

„ p 1 024 points that result. 
[0055] Because the GPS signal is extremely weak, the peak is unhkely to be detectable if 

Ifl processing integrates over only 1 millisecond. Consequently, results are improved if several 
; J milliseconds of input GPS signal are digitized, stored and processed. For example, using 4 or 16 

y 

' " milliseconds of data improves results due to gains that are achieved by integrating over these 
time periods to overcome the effects of noise. However many milliseconds are used, preferably 
the number is an integer multiple of 1 millisecond so that a complete Gold code sequence is 
observed in each millisecond of the signal collected. The Gold code (whose repetition interval is 
1 millisecond) should occur in the same place in each millisecond sample of the signal collected. 
[0056] Figure 6 is a schematic diagram of an embodiment of a system 601 for processing a 

multiple millisecond portion of the input GPS signal according to a third embodiment of the 
present invention. In Figure 6, a 4 ms portion of the input signal 602 is collected for processing. 
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As described above, the signal so collected can be any length, but is preferably an integer 
multiple of 1 millisecond. Input signal 602 is divided into 4 segments 604a5 604b, 604c and 
604d. 

[0057] An FFT processor 606 performs an FFT on each of the 4 segments. As described above, 

only the first 1,024 points are kept for further processing. The 1,024-point segments are 
segments 608a, 608b, 608c and 608d corresponding to input signal segments 604a, 604b, 604c 
and 604d, respectively. As would be known to those skilled in the art, this is a circular 
convolution, modulo 1 ms. Though circular convolution is sufficient in most cases, some 
processing gain is obtained by using a linear convolution where a transition occurs in the 
collected signal. This is described in more detail below. 
[005S| FFT segments 608a, 608b, 608c and 608d are multiplied by a fi*equency domain 

^ p representation of the Gold code 609 associated with a GPS satellite contributing to digitized 

input signal 602. As described above, fi'equency domain Gold code 609 is preferably pre- 
in computed. Alternatively, an FFT of the time domain Gold code is performed to obtain fi:*equency 
^0 domain Gold code 609. The multiplication proceeds point-by-point in multipUers 610a, 610b, 
610c, and 610d for each FFT segment and the product for each segment is stored in product 
612a, 612b, 612c and 612d, respectively. As described above, this operation is equivalent to 
time domain convolution. 

[0059] An inverse FFT algorithm 614 performs an inverse FFT on each product 612a, 612b, 

612c and 61 2d to yield convolutions 616a, 616b, 616c and 616d, respectively. Convolutions 
616a, 616b, 616c and 616d correspond to convolutions of the Gold code with input signal 
segments 604a, 604b, 604c and 604d, respectively. 
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If the carrier frequency is substantially correct, there should be a peak in each of these 
convolutions corresponding to where the Gold code is in each segment. It should be in the same 
place in each segment. Consequently, integration of the convolutions on a point-by-point basis 
should improve the ability to distinguish the peak from noise. A magnitude-and-sum algorithm 
618 is used to perform the integration. 

Convolutions 616a, 616b, 616c and 616d are applied to a magnitude and sum algorithm 
618. Magnitude and sum algorithm 618 performs a point-by-point magnitude and summing 
operation for each of the 1,024 points convolutions 616a, 616b, 616c and 61 6d. The sum of the 
points corresponding to the peak location should be greater than for other points which should be 
essentially noise values. A peak detector 620 determines the location of the peak by comparing 
each of the 1 ,024 points in the sum to a pre-determined peak threshold. If a value is greater than 
the peak threshold, the location where that value occurs is assumed to be the location of the peak. 
The peak location is the estimate of the code phase associated with the GPS satellite currently 
being processed. 

The correct carrier frequency is also determined as described above. If no peak is 
apparent, the frequency is shifted by the resolution of an FFT bin as described above. The half- 
bin analysis described above can be performed when processing multi-millisecond portions of 
the signal as well. In this instance, the half-bin processing is performed for each of the four 1- 
millisecond segments. That is, a complex exponential having a frequency equal to half a bin is 
multiplied by each of the 4-ms input signal segments 604a, 604b, 604c and 604d to shift the 
frequency by half a bin. Processing is repeated as described above, using the frequency-shift 4- 
ms input signal. In addition, if desired, zeros can be appended to products 612a, 612b, 612c and 
61 2d prior to applying inverse FFT algorithm to provide better bin resolution for peak detection. 
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Alternatively, as described above, frequency shifting can be accomplished by shifting 
frequencies in the frequency domain and shifting the Gold code FFT rather than the input GPS 
signal FFT. 

This processing is repeated for each GPS satellite contributing the composite input GPS 
signal. For each satellite, that satellites unique Gold code for each satellite. 

The foregoing processing assumed that there was no transition in the Gold code sequence 
of the collected portion of the GPS input signal. As described above, every 20 milliseconds the 
Gold code potentially inverts. The location of the transition in the 1 millisecond segment 
significantly effects the convolution of the Gold code for that particular segment. If the 
transition occurs at the beginning of the segment, there is a substantially negative peak. This 
does not affect the result since magnitudes are used in the summation. As the transition moves 
fiirther from the beginning of the segment, the peak degrades. It becomes negligible when the 
transition occurs in the middle of the segment. As the transition moves away from the middle of 
the segment, the peak that should result from the convolution increases in value. It returns to its 
maximum positive value where the transition occurs at the end of the segment. The worst case, 
therefore, is that the transition completely eliminates any contribution to the peak from the 1-ms 
segment in which it occurs. Consequently, if the transition occurs in the middle of a segment, 
then only 3 segments contribute to the peak where a 4 ms of input signal are integrated . 
Similarly, only 15 segments contribute to peak where a 16 ms of input signal are integrated. 
Thus, the transition is not likely to significantly degrade processing. 

The foregoing problem with the transition can be avoided by performing a linear 
convolution instead of the piecewise circular convolutions described above. An efficient class of 
techniques for performing such linear convolution is known generally as sectioned convolution. 
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In sectioned convolution, overlapping segments of the entire portion of the input signal stored in 
memory are convolved with the Gold code FFT stored in memory. The FFT of each overlapping 
section is taken and multiplied by the Gold code FFT. The inverse FFT of the result is 
calculated, and the overlap-save or overlap-add technique is performed using the result. Two 
well-known sectioned convolution techniques, the overlap-save and overlap-add techniques, are 
described in Lawerence R. Rabiner & BERNARD GOLD, THEORY AND Application of Digital 
Signal Processing, § 2.25 ("Sectioned Convolutions") (1975), which is hereby incorporated by 
reference herein in its entirety. 



convolution is a 4 ms linear convolution of the Gold code with the input signal. This signal is 

iu 

split into four 1-ms segments. The magnitude of each segment is calculated and the magnitudes 

ih 

for all of the segments are summed in a point-by-point fashion as described above. Because the 

iU 

2 sectioned convolution uses the entire Gold code sequence for each segment across segment 
IB boundaries, the difficulty described above with the transition where circular convolution is used 
is not present. Rather, the peak will be inverted where there is a transition. As described above, 
because magnitudes are used in the summing process, the presence of an inversion is inmiaterial. 
[0067] The process described so far uses non-coherent averaging. Averaging and integration are 



used synonymously herein. The location of the peak can be further refined using coherent 
averaging. A two-step process is employed for coherent averaging. In the first step, the 
above-described multi-segment technique is employed to obtain a coarse estimate of the carrier 
frequency and code phase. In the second step, coherent signal processing techniques are used to 
refine the obtained carrier fi-equency and code phase estimates. 



[0066L 



For the example, using a 4 ms portion of the input signal, the result of the sectioned 
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[0068] The coherent processing technique uses convolutions 616a5 616b, 616c and 616d. It is 

desirable to average these together by summation and then take the magnitude and perform peak 
detection. However, in order to perform such an operation, the frequency must be known very 
precisely. This is because, as you move across time in the input signal, any frequency difference 
between the input signal and the estimated carrier frequency manifests itself as a vector rotation, 
Le,, phase offset, in the data. This phase offset degrades the magnitude of the sums that are 
calculated. Over a millisecond, this vector rotation is probably not significant. However, over 
several milliseconds, this vector rotation can be detrimental to the averaging process. 
[0069] To account for the vector rotation, the frequency of the rotation is determined. This is 

^0 done by analyzing the complex values at the locations of the peaks (determined above) in each of 
the 1 ms convolutions 616a, 616b, 616c and 616d. The frequency is determined using any of a 
number of well-known techniques. For example, a four point FFT can be performed to 
determine the frequency of the sine wave. Accuracy is improved by the zero-appending 
technique described above. For example, better accuracy is obtained by appending 4, 12, 28 or 
hQ any number of zeros. Preferably, the size of the FFT is a power of 2 for computational 

efficiency. The FFT averages out the noise that is present in each individual value at the peak. 
Other curve-fitting techniques for determining the frequency of the sine wave using the complex 
values at the locations of the peaks would be apparent to those skilled in the art. 
[0070] Using this frequency offset, the rotation of the vector is compensated for in each of 

convolutions 616a, 616b, 616c and 616d. This is done by block multiplying convolutions 616b, 
616c and 616d by a complex exponential having a value of the sine wave of the frequency 
determined above at the location of the peak. Then convolutions 616a, 616b, 616c and 61 6d are 
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averaged. In this case, summing of convolutions 6163, 616b, 616c and 616d is performed prior 
to calculating the magnitude. This is possible due to the coherent nature of the processing. 

Figure 7 illustrates coherent processing in this fashion described above. Exemplary 
portions of convolutions 616a, 616b, 616c and 616d where the coarse convolution indicates the 
Gold code best aligns with the input signal for each millisecond of the input signal. These are 
peaks 701a, 701b, 701c and 70 Id. The phase associated with each peak is illustrated in phase 
representations 702a, 702b, 702c and 702d respectively. If the carrier frequency of the input 
si^al and the frequency of the locally generated Gold code were exactly matched, each of phase 
representations 702a, 702b, 702c and 702d would be identical. 

However, as shown in Figure 7, there is a rotation in the phase as time progresses, f.e., 
from convolution segment to convolution segment. The rate of rotation of this phase 
corresponds to the frequency offset between the carrier frequency of the input signal and the 
frequency of the locally generated Gold code. This frequency offset is estimated by, for 
example, taking an FFT, as described above. Arbitrarily, the phase of vector 702a is chosen as 
the reference. Using the frequency offset calculated above, the phase at the time of the peak for 
each of convolutions 616b, 616c and 616d is determined. The difference between this phase and 
the reference is the amount of rotation to apply to the respective convolution product to enable 
coherent averaging. The phase is applied by complex block multiplying the convolution segment 
with a complex exponential having a phase equal to the phase difference, from the reference. 

Further efficiency can be obtained by only processing the peak and one or a few more 
points on either side of the peak, rather than the entire 1-ms sample. Using only a few points is 
significantly faster than applying the process to all 1024 points. In the ensuing discussion, only 5 
points (2 points on either side of the estimated peak location and the point of the estimated peak 
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location) are used. It should be noted, however, that any number of points surrounding, and 
including, the estimated peak location can be used. 
[0074] After applying the rotation to the 5 points, a peak detection is performed to more 

accurately determine the code phase. There are many techniques for determining the peak. For 
example, curve-fitting techniques such as a parabolic fit, can be used to determine where the 
peak is given the points. 

[0075] Alternatively, a table lookup method is employed. In this method the Gold code is 

convolved with a perfect input signal at varying offsets, preferably over a chip length. 
Preferably, the offsets are varied uniformly over the chip length. Though any number of points 
could be used, for example, only 5 points may be used for computational efficiency. Any 
number of offsets can be used. For example, if 64 offsets are used, then the first set of values 

1 T\ 

corresponds to perfect alignment of the Gold code and input signal. The second set of values in 

; 3 : 

[ the table corresponds to a misalignment of 1/64^ of a chip length between the Gold code and the 
fg input signal. The table is preferably pre-computed and stored in non- volatile memory, as 

2 I 

described above for the FFT of the Gold code. 

[007^1^ Figure 8 illustrates a schematic for using a table 802 that stores pre-computed timing 

offsets of the input signal relative to the Gold code. The rotated 5 points in the vicinity of the 
peak 801 are correlated with each entry ENT 1, ENT 2 through ENT N in table 802 using 
correlator 804. The maximum correlation is chosen as the correct value of the code phase by 
maximum value selector 806. It is used to adjust the coarse code phase determined above from 
the magnitude and sum of convolutions 616a, 616b, 616c and 616d. 

[0077] Because the detection of the coarse estimate of the code phase and frequency is more 

computationally demanding than the interpolation required to fine-tune the estimates, the present 
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invention can be further optimized. This is because less data (for example 4 or 8 milliseconds) 
should be used in the initial step of determine coarse frequency and code phase due to the longer 
processing time requirements, while more data (for example 16 milliseconds) can be used in 
fine-tuning the frequency and code phase due to the shorter processing time requirements. 
However, it would be possible to apply the foregoing techniques to any amount of data. 
[0078] A variation of the present invention can be used where the SNR of the received signal is 

very low, for example, where the GPS receiver is in a building or other structure. In this case, 
the SNR may be too low to perform the non-coherent averaging technique to determine coarse 
frequency and code phase. What is required in this case is the ability to do coherent averaging to 
determine code phase and frequency. 
[0079"^ In this variation, processing is as described above to produce convolutions 616a, 616b, 

■p 616c and 616d. Then the coherent processing described above with respect to Figure 7 is 
I ^ perform on each point in convolutions 6 1 6a, 6 1 6b, 6 1 6c and 6 1 6d. Thus, rather than just 
In performing the coherent processing once at the coarse peak location, the coherent processing is 
=Q performed 1024 times (or however large the convolution size is). The resulting values are then 
i"^ compared to determine where the peak is. Once a peak location is determined, the peak location 
is refined as described above using curve fitting techniques or lookup tables. 
[0080] The foregoing disclosure of the preferred embodiments of the present invention has been 

presented for purposes of illustration and description. It is not intended to be exhaustive or to 
limit the invention to the precise forms disclosed. Many variations and modifications of the 
embodiments described herein will be obvious to one of ordinary skill in the art in light of the 
above disclosure. The scope of the invention is to be defined only by the claims appended 
hereto, and by their equivalents. 
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Further, in describing representative embodiments of the present invention, the 
specification may have presented the method and/or process of the present invention as a 
particular sequence of steps. However, to the extent that the method or process does not rely on 
the particular order of steps set forth herein, the method or process should not be limited to the 
particular sequence of steps described. As one of ordinary skill in the art would appreciate, other 
sequences of steps may be possible. Therefore, the particular order of the steps set forth in the 
specification should not be construed as limitations on the claims. In addition, the claims 
directed to the method and/or process of the present invention should not be limited to the 
performance of their steps in the order written, and one skilled in the art can readily appreciate 
that the sequences may be varied and still remain within the spirit and scope of the present 
invention. 
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